Spatial constraints underlying the retinal mosaics of 
two types of horizontal cells in cat and macaque 



Stephen J. Eglen 1,+ , James C. T. Wong 1 
Revision: 1.26 August 7, 2012 

i ■ 

q ■ Cambridge Computational Biology Institute 
CN ! Department for Applied Mathematics and Theoretical Physics 
bJ0| University of Cambridge 

Wilberforce Road, Cambridge CB3 OWA, UK 

. t Corresponding author. 

Phone: +44 (0) 1223 765761 
\ Email: S.J.Eglen@damtp.cam.ac.uk 

Postprint of Visual Neuroscience (2008) 25:209-214. doi:10.1017/S0952523808080176 



> 



Abstract 

Most types of retinal neurons are spatially positioned in non-random patterns, termed reti- 
nal mosaics. Several developmental mechanisms are thought to be important in the forma- 
tion of these mosaics. Most evidence to date suggests that homotypic constraints within a 
type of neuron are dominant, and that heterotypic interactions between diffe rent types of 
neuro n are rare. In an analysis of macaque HI and H2 horizontal cell mosaics, Wassle et al. 
(|2000l) suggested that the high regularity index of the combined HI and H2 mosaic might be 
caused by heterotypic interactions during development. Here we use computer modelling 
to suggest that the high regularity index of the combined HI and H2 mosaic is a by-product 
of the basic constraint that two neurons cannot occupy the same space. The spatial arrange- 



x 



ment of type A and type B horizontal cells in cat retina also follow this same principle. 
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horizontal cells, retinal mosaics, minimal distance model. 
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Introduction 



A defining feature for a type of retinal neuron is whethe r all neurons of that type tile the 
retina in non-random patterns, termed "retinal mosaics" (jCookl, 1998). This definition can 
also help us, together with other anatomical and physiological properties, determine whether 
a group of neurons should be classified as one type, or subdivided into several types. For 
example, cat beta retinal ganglion cells (RGCs) are classed into two types, the on-centre 
beta RGCs and the off-centre beta RGCs, partly because the mosaic of either the on- or off- 
centre neurons independently tiles the retina and each mosaic is much more regular than the 
combined mosaic of all beta RGCs ( Wassle et allll98lh . Furthermore, both cross-correlation 
analysis and modelling suggest that these two types of neuron are independent of each other 
in respect of positioning , as well as physio logical function, an d, hence, may dev elop inde- 
pendently dWassle et all Il98ll: lEglen etall l2005h . By contrast, IWassle et all (|2000h reported 
that for another pair of neuronal types, the HI and H2 horizontal cells in macaque: 

One would expect the nearest-neighbor distance of the combined mosaic to be 
smaller than that of the individual mosaics. The regularity index [defined in 
Methods, below], however, is comparable, suggesting that the HI and H2 cells 
are not arrayed completely independently. It is possible, that some interaction 
between their mosaics during retinal development creates this overall regularity. 
(|Wassle et alll2000L p597) 



In this report we use computer modelling to investigate whether the high regularity in- 
dex of the combined mosaic of HI and H2 neurons is a product of type-specific interactions 
between the two types, or whether it can be accounted for simply by anatomical constraints 
resulting from the two cell types occupying the same layer. To generalise this question 
slightly, and to evaluate more experim ental data, we wil l compare the spatial patterning 



of horizontal cells in macaque with cat (Wassle et ait [197; 



Methods 



Data sets Three horizontal cell fiel ds were analysed: fields A and B are from macaque 
(A: unpublished d ata; B: Figure 7 of Wassle et al. ( 2000h ): field C is from cat (Figure 12 of 
Wassle et al. ( 1978h ). To keep our notation concise (rather than claiming any equivalence of 
neuronal types across species), we denote type B cat horizontal cells as "type 1", and type 
A horizontal cells as "type 2" , in line wi t h pre viously-noted similarities of primate HI and 



other mammalian B cells ( de Lima et all 2005h . Fields were digitised, and the cell location 



taken to be the centre of each soma. Figure [T] shows an example real field along with a 
matching simulation, defined next. 



Bivariate d m i n model We have generalised the d m i n model (Galli -Resta et al., 1997) to sim- 
ulate the positioning of two neuronal populations within one field. Each type of neuron has 
its own homotypic exclusion zone (d\ or dj), but furthermore there is a heterotypic exclu- 
sion zone (d\2) to potentially allow for exclusions between the two types of neuron. (The 
subscript 12 refers to an interaction between two types of neuron, whereas the subscript 
i_i_2 used below refers to all neurons irrespective of type.) First, we count the number of 
type 1 and type 2 neurons (n\ and nj), and simulate an area A of the same size as the real 
field. To initialise the simulation, we randomly position n\ type 1 neurons and type 2 
neurons within A. Neurons are then repositioned randomly within the field subject to two 
constraints: that the nearest neighbour of the same type is greater than some distance {d\ 
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Figure 1: Real and simulated horizontal cell mosaics. Left: real mosaic (field A). Right: 
example simulation. Open circles denote HI cells, filled circles denote H2 cells; cells drawn 
assuming 10 pirn diameter. Scale bar: 100 }im. The simulated mosaic shows a close pair 
of H2 cells (halfway across, two-thirds up); such close pairs are rare but can occur when a 
homotypic exclusion zone is small. 



for type 1 neurons, &i for type 2 neurons), and that the nearest neighbour of the opposite 
type is greater than some distance d\j_. Each of the distances d\,di,d\i is a random variable 
drawn from a Normal distribution with a given mean and standard deviation. Random val- 
ues lower than a lower limit (5 ^m) are di scarded , to prevent implausibly small or negative 
dmin values. This birth and death process (|Ripley[|l977l;lEglen et all 120051) is repeated many 
times until convergence (typically after each neuron has been moved ten times). 



Null hypothesis Somata of both types of horizontal cell occupy the same stratum of the 
inner nuclear layer (INL). (In this study, we ignore the s mall population of displaced hori- 



zon tal cells that may be present in the ganglion cell layer (Sil veira et allll989l:IWassle et al. . 



bOQOh .) Our null hypothesis states that the developmental interactions between the two types 
of neuron that influence their positioning are limited to preventing somal overlap: any two 
neurons, regardless of type, cannot come closer than some minimal distance. In the context 
of our simulations, this implies that the range of d\j_ should match the range of typical somal 
diameters of the two types of neuron. 



Parameter estimates To fit one field, the free parameters in the model are the mean and 
standard deviation of the three exclusion zones. Th e homotypic exclusion zones {d\, dj) were 
estimated first by fitting a univariate d m i n model (|Galli-Resta et all Il997h separately to the 
type 1 and type 2 neurons. The size of the heterotypic exclusion zone, d\i, was assumed to 
be of the same order as the soma diameter of the horizontal cells, around 10 }im. Parameters 
were fitted by systematic searching over a range of plausible values. 



Assessing goodness of fit Two measures were used to quantitatively compare our model 
against the real data, the regularity index (RI) and the K function. The RI is computed by 



3 



field 


ni 


n 2 


width x height (jim 1 ) 


d\{}lTS\) 




rfi2(^m) 


A 


187 


82 


400 x 402 


22 ± 4 


40 ± 10 


11 ± 3.0 


B 


206 


86 


298 x 300 


21 ± 4 


32 ± 8 


12 ± 2.5 


C 


300 


85 


723 x 1194 


65 ± 12 


72 ± 8 


14 ± 3.0 



Table 1: Parameters for each simulation. The number of type 1 and type 2 neurons {n\, 
n-i) and the field size matches the values from the real field. The diameter of each exclu- 
sion zone {d\, &i> dyi) is drawn from a Normal distribution, listed here as mean ± standard 
deviation. 



measuring the distance of each non-border neuron to its nearest-neigh bour, and then di - 
viding the mean of this distribution by its standard deviation (|Wassle & Riemannl, 119781) . 
(Nearest-neighbour distances of neurons at the border of a field are excluded as those dis- 
tances are unreliable. A neuron is excluded if its Voronoi polygon touches the boundary of 
the field. This exclusion criterion accounts for the small differences in RI between our work 
and those previously reported.) We measure three RI values: RIi (distance of each type 1 
neuron to nearest type 1 neuron), Rfy (like Rl\, but for type 2 neurons), and RI\ + 2'- (distance 
of each neuron to nearest other neuron, irrespective of type). 

For one population o f neurons, Kit) measures the n umber of cell pairs within a given 
distance t of each other (jRipleyi, [19761: lEglen et al. . 2005h . For plotting purposes, we show 
L(f) = [K(t) I n] 1 ^ 2 . This transformation discriminates between exclusion (L(f) < t), clus- 
tering (L(f) > t) and complete spatial randomness (L(f) = t). We measure four L functions: 
L\\ pairs of type 1 neurons; L^. pairs of type 2 neurons; Li + 2'. pairs of neurons of either 
type. Finally, L\2 measures the cross-correlation, by constraining cell pairs such that one 
cell is type 1 and the o ther is type 2. Full details of these measures are given elsewhere 
(|Eglenetal.l.E003l.l2005h . 

To quantitatively evaluate the goodness of fit of the model to the real data, each sim- 
ulation was run 99 times with the same parameters, but from different initial conditions. 
Informally, if the measure from the real data falls within the distribution of observed values 
from the simulations, then the model fits the data. This can be quantified with a p value 
using a Monte Carlo ranking test. A test statistic (Tj) is measured for the K function of the 
real mosaic (z = 1) and for each simulated mosaic (i = 2 . . . 100). A p value is then calculated 
by dividing the rank (smallest first) of T\ by 100. P values greater than 0.95 indicate a signif- 
icant dif ference at the 5% level between model and data. Full details of the test statistic are 
given in lEglen et all (|2005h . 

Computational modelling and analysis wa s performed in the R environment, using the 



splancs package and Voronoi domain software (R Development Core TeamLl2007i;lRowlingson & Diggle 
1993: lFortune].ll987h . as well as custom-written routines. The code is available from the au- 



thors upon request. 



Results and Discussion 

Table Q] lists the parameters used for each bivariate dmin simulation. The homotypic exclu- 
sion zones (di, ^2) were independently fitted to each mosaic, whereas the heterotypic exclu- 
sion zone (di2) was set to just prevent neurons of opposite type from occupying the same 
space in the inner nuclear layer. (Mean values of d\2 reported in Table [T] are slightly higher 
than estimates of somal diameter, suggesting that both cell bodies and some initial portion 
of the primary dendrites contributed to steric hindrance between neuronal types.) Figure [2] 
shows that for each field, the model generates mosaics that quantitatively match the real 
mosaics, as assessed by both the RI and the L functions. In three (out of twelve) cases the 
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goodness of fit p value is greater than 0.95, indicating that formally there is a significant dif- 
ference between data and model. Two of these cases concern both type 1 and type 2 neurons 
from cat retina (field C). Discrepancies in these two cases are apparent over small distances 
(less than 10 pirn between neurons of opposite type), and may be due simply to difficulties 
in reconstructing the position of pairs of opposit e-type neurons tha t seem to overlap in the 
field from the original publication (Figure 12 of IWassle et all 119781) . Small errors in deter- 
mining neuronal position are likely when considering the relative size of individual neurons 
with the size of the sample field. Overall, however, the L functions for the data fit within the 
confidence intervals of the model, suggesting that any disagreements between model and 
data are quite small. 

Regularity of a d m i n mosaic is influenced by both neuronal density and the distribution of 
exclusion zone diameters. For both macaque fields, the median RI is higher for the simulated 
type 1 mosaics than for the simulated type 2 mosaics; the opposite is true for the cat field 
(p < 0.001 in each of three cases, Wilcoxon rank sum test). The high RI of simulated cat 
type 2 mosaics, matching the observed data, is due to the relatively low s.d. in the type 
2 exclusion zone (Table [T]); if this is increased (e.g. from 8 pm to 16 f/m), the median RI2 
decreases to 4.5 (data not shown), below that of the type 1 mosaics. 

The RI of the combined type 1 and 2 mosaic (RIi + 2) in each of our fields is typically 3-5, 
matching the values observed experimentally. At first glance, this might seem quite high, 
especially compared to a theoretical e xpecte d value of around 1.9 for cells (of infinitesimally 
small size) arranged randomly (|Cookl,ll9 96). Our model tells us that this high RI is simply a 
by-product of superimposing two regular, but independent, mosaics with somal exclusion. 
This conclusion can be supported in two ways. First, by setting d\i to zero, we eliminate 
all heterotypic interactions. (This allows for neurons of opposite type to become arbitrarily 
close to one another, which is of course not realistic, but allows us to specifically test the 
impact of removing all heterotypic interactions.) Figure [3]A. shows that RI1+2 drops con- 
siderably to a median of around 2.5, whereas the fits to Rl\ and RI2 remain good. In this 
case, since there is no positional constraint between opposite-type neurons, the L function 
for the random simulations follows the theoretical expectation L\i{t) = t. The deviation 
between real data and simulations is apparent up to at least 20 jiva (Figure |3p), as observed 
by the L function for the real curve dropping well below the confidence intervals from the 
simulations. 

The second line of evidence to explain the high RI of the combined (type 1 and type 2) 
mosaic is suggested by examining the fraction, /, of the retinal area occupied by the cell 
bodies. This can be estimated by / = {{n\ + n2)7ir 2 ) /\A\ where r = 5 jim is an estimate of 
radius of a horizontal cell soma, n\ and nj_ are the number of type 1 and type 2 horizontal 
cells, and | A | is the area of the field. Table |2] shows that the fraction of occupancy (/) corre- 
lates with the regularity of the combined mosaic {Rl\ + 2)- Furthermore, Table |2] also shows 
that the mean number of trial cell positions rejected due to infringement of the heterotypic 
constraint also correlates with the RI, even when normalised for the number of potential 
pairwise heterotypic interactions. In this light, the cat and macaque mosaics are generated 
by the same mechanism, and the lower regularity of the combined mosaic in cat is due to 
the smaller effect of somal exclusion. 

Our bivariate d m i n model is open to criticisms of biological plausibility. As previously 
noted, exclusion models show us that local interactions are sufficient to g enerate regular 
patte rns, but do not inform us on how these interactions are mediated (|Galli -Resta et al., 



19971 ). Here we would suggest that the homotypic exclusion zones are mediated by horizon - 



tal cell processes, perhaps driving lateral migration during development (Re ese et al. . 1999). 
The heterotypic interactions however are simply the result of steric hindrance between cell 
bodies and primary dendrites. 
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Figure 2: (Colour online.) Goodness of fit between data and model. Each column compares 
one field (A, B: macaque; C: cat) with simulations. In the regularity index (RI) plots, the 
thick green line indicates the RI of the real data for either type 1, type 2, or all (1+2) neurons. 
Black dots are RI values from 99 simulations, together with their median (thin black line). 
For each of four L function plots for a field, the L function for the real mosaic is shown as 
a solid green line and dashed black lines indicate 95% confidence intervals of simulations. 
The p value is the goodness of fit between model and data. 
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Figure 3: (Colour online.) Results of bivariate d m in simulation for field A with d\2 = 0, and 
other parameters as listed in Table [TJ Results are presented as in Figure |2] 



field 


RI 1+2 


f 


rejects 


Wpairs = n 1 xn 2 


rejects / n pairs 


C 


3.4 


0.04 


242 ± 10 


25500 


0.01 


A 


3.9 


0.13 


940 ± 32 


15334 


0.06 


B 


4.9 


0.26 


7198 ± 282 


17716 


0.41 



Table 2: Estimates of the fraction / of sample retinal area occupied by all horizontal cell 
bodies in each field and incidence of heterotypic constraint enforcement. Rows are sorted in 
order of increasing regularity index of the combined mosaic. The mean (± s.d.) number of 
times (per sweep) that the heterotypic constraint was broken, rejects, was counted over 99 
simulations. The final column shows the mean number of rejects divided by the number of 
pairs of opposite type neurons. 
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In conclusion, results from our computational model suggest that the high RIs observed 
for combined HI and H2 mosaics in macaque are simply a by-product of the two mosaics 
being positioned in the same stratum of the INL; the mosaics may be developmentally in- 
dependent in all other respects. This result agrees with our earlier work on beta RGCs 
(Egl en et"aD, 2005), as well as other studies demonstrating a lack of spatial c orrelations be- 
tween many pairs of retinal neuronal types (IMackl, 120071: iR ockhill et all 120001) . Exceptions to 
this finding are rare (|Kouyama & Marshal. Il99ft I Ahnelt et allbOOCh . 
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